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We show through a variational calculation that in a large 
range of parameters the paramagnetic to ferromagnetic transi- 
tion in colossal magnetoresistance compounds is accompanied 
by a collapse of polaronic lattice disorder in addition to that 
of spin disorder. The spin-lattice disordered state is shown 
to be localized and the ordered state itinerant as observed. 
The observed dramatic change in the diffuse scattering at the 
transition as well as the isotope effect are also explained. 

PACS numbers: 72.15.Gd, 75.30. Kz, 75.30.Vn 

While the basic physics of the colossal magnetoresis- 
tance (CMR) phenomena follows |||-^ from the 
Double- Exchange (DE) model, several anomalous 
results have emerged in a number of recent experi- 
ments ]lO|-p^, which require additional considerations: 
(i) the high-temperature paramagnetic phase in these 
compounds is often characterized by insulating behav- 
ior; (ii) the ferromagnetic phase transition observed in 
these materials is often first order; (iii) the paramag- 
netic phase is characterized by an intense diffuse scat- 
tering, which disappears with temperature in a manner 
correlated with the spontaneous magnetization; and (iv) 
the transition temperature displays huge isotope effects. 
These phenomena are related in that when there is lit- 
tle diffuse scattering, the transition is second order or 
weakly first order, the resistivity above the transition is 
metallic and there is only a small isotope effect. 

The diffuse scattering experiments reveal inhomoge- 
neous lattice configurations in the paramagnetic phase 
which anneal dramatically in the vicinity of the CMR 
transition ||l0| , [rT[| . Spin and lattice "disorder" in CMR 
compounds are thus closely related. We believe this pro- 
vides a important clue to the essential physics of the 
colossal magnetoresistance. In the present Letter, we 
present a systematic method to describe this effect and 
its consequences. 

The problem we solve also has a more general context. 
There are two different ways in which electron-lattice in- 
teractions are considered. One, which we will call the 
polaron or the Holstein approach, first deals with 
the problem of a localized electron interacting with the 
lattice to form a polaron, then the hopping from site 
to site is introduced, and the polaron Bloch-states are 
formed. ||lj,|f^ In the other approach, which may be as- 
sociated with Migdal, one first forms the Bloch-states 
of the electrons, and then treats the electron-phonon in- 



teractions perturbatively in Q/Ep, where 9 is the De- 
bye temperature and Ep is the Fermi energy. These two 
methods yield vastly different results, the electron self- 
energies being typically of the order of the Fermi energy 
in the first case, and O (8) in the second. The transi- 
tion behavior between these two limits is however not 
yet adequately understood, and has recently attracted 
considerable attention. |p7|] 

The model we consider is defined by the Hamiltonian 
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The first three term represent the standard double- 
exchange (DE) model with transfer integrals = 

i°.cos^ 



which depend on 6ij, the mutual orientation 



of spins at sites i and j. The remaining terms in Eq. (|1|) 
describe the electron-lattice interaction, and a model lat- 
tice Hamiltonian. Due to the large difference between the 
ionic radius of Mn^^ and Mn**^, the breathing mode of 
the Oe octahedra is expected to be the most important, 
although asymmetric or Jahn- Teller modes may also con- 
tribute. ]l8[ | The dynamics of the oxygen octahedra sur- 
rounding manganese ions, is represented in our model by 
localized classical oscillators with a scalar displacement 
Ui from some reference point and the corresponding mo- 
mentum Pi. The principal physics underlying the CMR 
phenomenon, can be uncovered using one coordinate per 
each Mn ion. Similar models for the CMR problem have 
been considered before. [0 

The essential physical point we focus on is the com- 
petition between the electrostatic energy, which favors 
large lattice distortions provided the electrons are local- 
ized, and the gain in kinetic energy due to itinerancy of 
electrons which suppresses lattice distortions. This com- 
petition is self-consistently reinforced by spin-disorder fa- 
vored by entropy and spin-order favored by the kinetic 
energy. (Extrinsic disorder tilts this competition in favor 
of localized states.) Accordingly we consider a variational 
density matrix: 
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where /i is the chemical potential. The energies Sk are 
determined variationally, and 
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where the operator |j, u, i9) (j, u, t?| represents a projection 
on a state with the lattice distortion u and spin direction 
'd at site j, while ui and U2 are the lattice distortions at 
"occupied" and "empty" sites respectively. The proba- 
bility distribution P is also determined variationally as 
described below; different limits of it describe param- 
agnetic localized electronic states with distorted lattice 
configurations or itinerant ferromagnetic states with uni- 
form lattice configurations. 

Neglecting local spin-lattice correlations we represent 
P as : 

P{ui,U2;du'd2) ^Pu{ui,U2)PAdi)P^{d2) (4) 

Substituting (jj) into the free energy obtained from Eq. 

and taking the functional derivative with respect to 
the lattice distortion distribution P„, after a straightfor- 
ward but lengthy calculation we find that the energies Sk 
correspond to the eigenstates of the Hamiltonian 



F. = (cos-)(exp(^-^ 



where Su = U2 — ui. 

The variational procedure yields the physically trans- 
parent result for the lattice distortions distribution 



Pu{ui,U2) (X exp{-E[ui,U2] /T) 
so that the effective free energy 

E[ui,U2] 
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F[p^] =-riog 
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where 



E[ui,U2] = / dep (e) e — Xuix + {ui,U2) 
and the lattice energy 



El {ui,U2) = - {xuj + (1 - x) ul) 
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X is the number of electrons per site, p (s) is the electronic 
density of states (DOS). The first term in Eq. (||), 

dep (e) e ~ (1 — a;) W, 

-00 ^ 

where cq is a constant which depends on the functional 
form of the density of states (e.g. for rectangular DOS 



Co = 1, while for the Gaussian functional form of p (e) we 
find Co ~ 2.26 and for "elliptical" DOS co w 1.7). The 
bandwidth W depends both on the spin distribution and 
the lattice distortions: 



W :^Wg exp 
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where We = Wo (cos |). For large Su the Eqn. ( [lO| ) 
represents the suppression of the bandwidth due to small 
polaron formation. 

The critical temperature Tc for the CMR transition is 
typically much smaller than the electronic cohesive en- 
ergy. Therefore, to adequately describe phenomena in 
the vicinity of the transition and below it we need con- 
sider only the limit T <C E (ui,U2). Then the integral 
over the distortions ui and U2 in the effective free energy 
Eq. (|^) can be calculated by the method of steepest de- 
scent. Changing the variables to u = i (iti -|- U2) and Su, 
we find 
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where Eu depends only on Su : 
Co 



Eu (Su) = X {I — x) 



-W (Su) - XSu + -Su'^ 
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(12) 



and the index j labels different minima of the effective 
energy Eu- 

Depending on the ratios of the electron bandwidth Wq 
to the polaron energy Ep = A^/fc and the phonon energy 
hujo, the function Eu (Su) has either one or two minima. 
The corresponding phase diagram is presented in Fig. 
|l|. In CMR compounds, the bandwidth Wq ^ hioo, and 



therefore the summation in Eq. (11) always includes two 



and only two terms. The minimum at small distortions. 
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which we will refer to as the "Migdal 



kcoWo (cos I) ■ 

minimum" , corresponds to extended electron eigenstates 
with renormalized effective mass with the energy 
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The other minimum, which we will refer to as the "Hol- 
stein minimum", occurs at m*-^^ ~ X/Mujq, and corre- 
sponds to small polarons with the energy 
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Substituting (|l|) and (|l|) into Eq. (0), we obtain 



(14) 
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FIG. 1. The "phase diagram" corresponding to the en- 
ergy Eu (Su), which represents the number of energy minima 
as a function of Wo/hujo and Wo/ (A^/fc). The green area 
corresponds to a singe minimum of Eu, while the blue and 
the red regions represent the parameters when the energy has 
two minima with the lowest being either the Holstein mini- 
mum (the blue region) or Migdal minimum (the red region). 
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where -7 ~ W 1 4- 2cox{i-x) Wg_ 
wneie 7 _ y i -t- i^^^f^^^ j^^^ 

The equations (|l^), ( |T^ ) and ([l^ ) define a variational 
problem for the spin distribution function , which can 
be solved by a straightforward minimization of the free 
energy. The corresponding procedure is similar to the 
one used for the standard DE model and described in 
detail in Ref. |l^. We generally find a phase transition 
between low lattice distortion ferromagnetic phase at low 
temperatures and high lattice distortion paramagnetic 
phase at hight temperatures. For a strong electron-lattice 
couphng, fcWo ^ , the transition is of the first kind, 
and is accompanied by the virtual collapse of the lattice 
distortion 5u, as shown in Fig. |^. Note the remarkable 
similarity of the predicted behavior to the observations 
in recent experiments [ p^pd] , as illustrated by the inset. 

In the other limit, when the electron-lattice coupling 
is small ( i.e. kWhtoo » A^), our model corresponds to 
the standard double-exchange, and the transition is of 
the second kind. In fig. |^ (a) we plot the change of the 
lattice distortion Su = U2 — mi at the phase transition, 
Au\c, as a function of the model parameters, A^/fcWo and 
Wo/huJo, in false-color representation. There, the "deep 
blue sea" for small electron-lattice coupling corresponds 
to the second-order phase transition, where Au\c = by 
definition. 

The phase transition changes from the first order kind 
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FIG. 2. The magnetization M/Msat (dashed blue line) 
and the lattice distortion Au in units of A/fc (green line) as 
functions of temperature T (in units of Wo), for the carrier 
concentration x = 0.3. The inset shows the experimental 
temperature dependence of diffuse scattering intensity due to 
lattice distortion from Ref. [^. 



to the second order when the ratio of the polaron en- 
ergy Ep = A^/fc to the electron bandwidth Wq is near 
unity. This is in a remarkable agreement with experi- 
mental data, as these quantities in lanthanum mangan- 
ites are believed to be of the same order, and both first 
|l0| , p| and second-order phase transitions are indeed 
observed in experiments depending on the composition 
(which is likely to affect the electron-lattice coupling and 
the bandwidth). 

An increase of the electron-lattice coupling leads to 
suppression of the metallic phase, which manifests itself 
in smaller values of the critical temperature. This behav- 
ior is illustrated in the bottom panel of Fig. ^ , where 
we present the false-color representation of the depen- 
dence Tc as a function of both parameters of the model, 
X^/kWo and Wo/huja. 

An important feature of the proposed theory is the 
strong dependence of the critical temperature of the fer- 
romagnetic transition on the effective ionic mass M . This 
behavior is illustrated in Fig. |^ where we plot the criti- 
cal temperature (in units of the "bare" bandwidth Wq) 
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M. As follows from Fig. W, a 10% 



change in the oxygen ion mass which corresponds to the 
replacement of ^^O by ^^O, could account for a varia- 
tion of critical temperature by 20%, which explains the 
recently observed giant isotope effect in lanthanum man- 
ganites. |l^] 

Note the remarkable "anti-correlation" of the change 
of the lattice distortion at the transition with the critical 
temperature, as seen from Fig. |^. This finding is again 
in agreement with the experimental data : the 

transition of the second order kind is generally observed 
in the compounds with relatively high values of Tc, while 
a second-order transition is found in materials with a 
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FIG. 3. The change of lattice distortion at the transition 
in units of Up = X/k (top panel) and the critical temperature 
in units of To = — x)Wo (bottom panel), as functions 

of X^/{kWo) and Wo/hujo, in false-color representation. The 
"deep blue sea" in the top panel represents the parameters 
for the second-order transition when Aw = 0. The carrier 
concentration x = 0.3. 

smaller value of critical temperature. 

The dependence of Pu{ui, U2) on Am in our ansatz for 
the distribution function (Q) introduces the essential cor- 
relation of lattice distortion with electronic occupation. 
If this correlation is ignored by approximating Pu{ui, U2) 
in a mean-field fashion by p{ui)p{u2) , the important re- 
sults of our theory presented in Figs. ||-^ which include 
the collapse of the lattice disorder at first-order ferromag- 
netic transition and strong isotope effect in the critical 
temperature, cannot be obtained. 

We conclude with a discussion of the effect of intrin- 
sic disorder. The intrinsic disorder, which can be rep- 
resented in the model Hamiltonian (|^) by a term of the 
type ^iViclci, favors less ordered paramagnetic phase 
and therefore leads to a decrease of the critical temper- 
ature. This effect was considered in detail in Ref. p9| . 
More important is the effect of the intrinsic disorder on 
the transport properties of the CMR compounds. In 
high quality samples with relatively weak disorder, the 
"Migdal minimum" describes the extended states and is 
responsible for the metallic behavior below the transition 
temperature. 

If the electron-phonon coupling is substantial, the situ- 
ation is dramatically different above the transition where 
the dominant contribution is given by the "Holstein min- 
imum". There, the effective bandwidth is exponentially 
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FIG. 4. The critical temperature as a function of 
Wo/hwo ~ VM (for a fixed value of X^/{kWo) = 1.). 

suppressed due to lattice distortions 

Wh (X Wo exp f-^^vA (16) 
V n^ivo J 

where A^/fc ~ Wq 3> hujo (in manganites Wo/hujQ ^ 50) 
and V = k 6u/X is the dimensionless lattice distortion 
which in the paramagnetic phase is of the order unity 
(see Fig. ||). Therefore Wh ^ Wq, and small intrinsic 
disorder in the CMR compounds will lead to the local- 
ization of the Holstein polarons. On the other side of the 
ferromagnetic transition, small intrinsic disorder alone, 
because of the vanishing spin as well as polaron disorder, 
is insufficient to localize the electronic states. So metal- 
lic conduction occurs below the ferromagnetic transition, 
while an activated form of resistivity is obtained above 
the transition, in agreement with experimental data. 

In contrast, for low electron-phonon coupling con- 
stants, as evidenced for example by high (second-order) 
transition temperatures, in pure samples the resistivity is 
expected to be "metallic" on both sides of the transition. 
Experimental results are consistent with these observa- 
tions. 

To summarize, in this paper a systematic method has 
been found to isolate and calculate the principal effects of 
electron-lattice interactions in strongly correlated mate- 
rials of the double-exchange family. By including effects 
of the correlation of lattice vibrations with local elec- 
tron occupation, we have shown that for the same rea- 
son that the itinerant state sweeps away spin-disorder, 
it also sweeps away lattice disorder. This has been ac- 
complished by a variational form for the couple electron- 
lattice that extrapolates from the Holstein form to the 
Migdal form as the self-consistently determined band- 
width changes from a value small compared to polaron 
binding to a value large compared to it. The results of 
the calculations explain many of the peculiar features of 
the experimental results of CMR compounds. 
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